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Abstract 

The article is devoted to the dynamics of systems with an anomalous scaling 
near a critical point. The fractional stochastic equation of a Lanvevin type 
with the (p^ nonlinearity is considered. By analogy with the model A the field 
theoretic model is built, and its propagators are calculated. The nonlocality 
of the new action functional in the coordinate representation is caused by 
the involving of the fractional spatial derivative. It is proved that the new 
model is multiplicatively renormalizable, the Gell-Man-Low function in the 
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the framework of the e-expansion for a superdiffusion is established. 
Keywords: 

fractional processes, renormalization group, stochastic evolution equations, 
critical scaling, nonlocal lagrangian 



Email addresses: zlokor88@gmail.com (L. A. Batalov), AsjaStallOgmail.com (A. 
A. Batalova) 



Preprint submitted to Physica A: Statistical Mechanics and its ApplicationsMarch 4, 2013 



1. Introduction 



The models based on the stochastic processes with continuous time are 
widely used in modern theories of phase transitions, an atmospheric tur- 
bulence and electromagnetism , {2!] , • From the thermodynamic point of 
view these models are characterized by the existence of the temperature Tc 
such that observed experimental results for T — t- Tc — and T — )■ + differ 
significantly {4]. 

The object of our investigation are systems governed by stochastic equa- 
tions in the vicinity of the temperature Tc called the critical point. Problems 
of an equilibrium statistical physics and a thermodynamics without an ex- 
plicit time-dependence are called static. Obviously in such systems there are 
processes the description of energy-dissipation. The description of a critical 
behaviour in the dynamic models is based on the analogous static case, but 
it is more complex technically jsl. In critical dynamics only the models of a 
Langevin type are considered, i.e. the models with the first time derivative 
on a field. 

Among these models the Ginzburg-Landau dynamic model (the model 
A in the Galperin and Hoenberg classification scheme jsl) proposed by Lan- 
dau and Khalatnikov for description of a superfluidity is the most interesting 
. The critical dynamics of the model A was investigated using a dimen- 
sional regularization and RG- methods based on the e-expansion [7| . The RG 
method allows to prove the existence of the critical scaling (it is a power- 
law dependence near a critical point of dynamic correlators) and to obtain 
corrections for the main exponent jsj. 

For the isotropic system which consists of particles moving under the 
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influence of tlie random Gaussian distributed forces the base of theoretical 
analysis is the diffusion equation. 

For a classical diffusion in the d - dimensional space a particle travels 
for the time period t the distance At'^'^/^, where A is the Onsager kinetic 
coefficient, which here can be understood as the diffusion coefficient. This 
dependence is called a normal scaling. However, the power-law dependence 
of the displacement with the critical exponent more or less than 1/2 (an 
anomalous scaling) is observed in some physical systems jsj. 

For such FDE the corresponding scaling exponent is equal a/a, and frac- 
tional derivatives are the generalization of integer derivatives. It allows to 
apply for a solution of linear equations standard methods the Laplace and 
Fourier transforms. 

In the particular case a = 1, a = 2 we obtain the standard equations 
of the Brownian motion. An exact solution of the linear fractional diffusion 
equation was found for many processes. Some of them generalise the solutions 
of the standard diffusion equation with the exponential functions replaced 
by the Mittag-Leffler function lo|. Moreover the class of exact solutions 



of non-linear diffusion eq 
very much last time 
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nations and the Fokker-Plank equations expanded 



12|. 



The RG methods and the ideas about description anomalous scaling 
with fractional differential equations (FDE) with derivative orders a (a time 
derivative) and a (a spatial derivative) were developed in the same time. 

In our opinion the application of the renormalization group to the solu- 
tion fractional differential equations will produce good results. This sugges- 
tion bases on the observation, which shows us, that much of the systems 
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demonstrated the critical behaviour have the scahng functions with anoma- 
lous exponents. Governing equations of these systems are fractional (for 



example a developed turbulence 



13|). In our paper the technique of the RG 



method generalisation on dynamic models with fractional spatial and time 
derivatives is proposed. 

We have developed a technique of a renormalization for an arbitrary 
stochastic problem with the nonlinearity. The one-loop diagram for the 
most essential model of the critical dynamics - the model was calculated. 
We showed, that an analogy with statics took place only for a = 1 and the 
independence on a and m. 

This technique is convenient in applications, it is used in the calculation 
of higher orders, because it is_based on the well-developed quantum-field 
renormalization-group theory 



14|. However, it will be showed in the paper, 
that for the general analysis and determination of critical regimes the 1- 
loop approximation is enough. Before now the RG method was applied to 
a nonequilibrium system usually with an anomalous diffusion in the Wilson 
recursive equations form 5|, isl- We emphasize, that in our work it was used 
the quantum-field renormalization. The critical scaling will be proved, the 
main exponential function and their dependence on a and a will be found. 

The article construct as follows: in section 2 the stochastic task of the 
critical dynamics is manifested, and its formulation is given in the field the- 
oretic terms. The section 3 is devoted to the calculation of propagators of 
this model. In section 4 the model is analysed thoroughly and a transi- 
tion to the renormalized functional of action in the logarithmic dimension of 
space is carried out. Below we derived the RG-equations and calculated the 
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Gell-Mann-Low function in the one-loop approximation. 



2. Formulation of problem 

We consider a nonequilibrium system described of the scalar function 
(p{'x,t), where t is time, x - ci - dimensional radius vector. All calculations 
demonstrated in this section are easy generalized on the case, where field ip 
is multi-components. 

The critical dynamic of field <f{'x, t) is described by the stochastic Langevin 
equation: 



where V{(p) = 6S/6(f - variational derivative of the static action functional 
S{ip) of time-independence free field ipi^), A > - Onsager transport coeffi- 
cient, F - random external force. For this a Gaussian distribution with zero 
- mean (F(x)) = and covariance 



is assumed. 

The equality of the double Onsager coefficient and the correlator power 
provides an agreement with the static case described by the functional S{ip). 
It is consequence of the Fokker-Planck equation for the field distribution 

n 

function P{(p) [7|. 

The functional potential V{(f) includes a linear differential operator of 
the second order (Laplacian) for different generalizations of the Ginzburg- 
Landau functional (further - the ip^ model). Therefore the linear part of the 



(1) 



(F(x, t)F(x', t')) = 2A(5(x - x')5(t - t') 



(2) 
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stochastic equation (ITl) is a linear combination of the first derivative time 

R 

operator and the second derivative space operator [6]. 

The normal scaling is concluded from Eq. ([1]) in the leading order go = 
and V{ip) = A(p: (x^) = At^'^^, where H, called the Hurst exponent, which 
is equal 1/2. 

We can use Eq. ([1]) included instead of the partial time and spatial 
derivatives the fractional derivatives of orders a and a related to the Hurst 
exponent H = a/a for description of the anomalous diffusion processes with 
same parameter A. For representation of subdiffusion processes < H < 1/2 
the partial time derivative could be changed by the fractional derivative of 
order a (0 < a; < 1), and in the case of super diffusion 1/2 < H < 1 
laplacian could be changed by the differential operator of order a. In models 
of turbulent transport the Hurst exponent is greater than 1/2 and is close to 
3/4. 

Therefore we will find so generalization of the problem ([T]) , that allows to 
use the fractional derivatives. It is convenient to use the time derivative of 
order a '^Df in the Caputo representation with the lower terminal a, which 
is usually chosen from the condition of zero antiderivative [3]: 

V{n-a)J, (t-r)"+i-" 

where n is the nearest integer bigger than a. This representation provides 
the Laplace transform as following: 

m=0 

As a generalization of the Laplacian, it is convenient to define the operator 
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in momentum representation, in other words, we will work with its Fourier 

image. The word "momentum" here can be understood as a wave vector. 

If (T - order of the spatial derivative and (— A)'^/^ - the so-called frac- 

□ 

tional Laplacian acting on the function of the Schwartz class [17|] /, so that 
(_A)-/2/ = h, then 

^(k) = r/(k), k = |k| (5) 



If < cr < 2 we can also use the integral representation [17l |: 

(-A)'^/V(x) = C{n, a)P.V. j ^l^\~J^] dy (6) 

where P.V. stands for principal value and C(?t,, a) = 2'^~^crr(((i+cr)/2)/(7r^/^r(l- 
(7/2)) is a normalization constant. Space R'^ has standard euclidean norm 

This allows to apply to diffusion equation Laplace or Fourier transforms. 



The diffusion equation arising in this way is exactly solvable model 18|. 

^AV(x,t) = -A(-A)'^/V(x,t) (7) 
supplemented by the following initial and boundary conditions: 

¥?(x,0) = 5(x), lim (/?(x,t) = 0, t > (8) 

|x|— >oo 

where A > - the diffusion coefficient, 0<a<l, 0<cr<2, v^(x, t) - scalar 
field, qD^ - the defined on positive half-line Caputo derivative of order a. 
The Green function of Eq. (171) with the initial and boundary conditions 

n n 

dH]) in the time - momentum representation is equal [ld],[18| Ea{—\k°'t°') 
where Ea{z) - the Mittag-LefHer function, which is defined usually by the 
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following series: 

n=0 ^ ' 

In the case of real and positive a the series converge for all values of the 
argument 2;, thus the Mittag-Leffler function is an entire function. 

There is the simple and universal technique to transform the stochastic 
equation dl]), (E]), (IHl) to a field theory by doubling the number of fields 



191]. It yields an explicit form of action functional >S'($), where $ = (v?, V'') 
is a new two-component field. We plan to generalise this technique on the 
fractional derivative problem. Note, that the translational invariant pertur- 
bation theory with zero condition for y9, if t — )■ —00 is considered in critical 
dynamics, than time integration is performed on the entire axis in all terms. 
The initial condition ([8]) means, that integration is performed on the values 
t > only. 

We will consider the following generalization of Eq. ([1]) 

SAV = >^V{ip) + F (10) 
with a non-linear potential of the nonlinearity degree m > 0: 

Viif) = -i-Ay/^cp - Tocp - g,v^/m\ (11) 
with the initial and boundary conditions 

lim (^(x,t)=0, ^(x, -oo)=0 (12) 

|x|— >-oo 

and the conditions on orders of derivatives: 

0<a<l, 0<(j<2 (13) 
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and we assume the positivity of temperature and coupling constants: tq > 
0, go > 0. 

The power of correlation function (|2]) remains the same and it is equal 2A, 
because in the proof of static and dynamic theory consistency the condition 



V{(p) = 6S/6ip is used It will be proved later. In the case of the Cauchy 
problem ([8]) the equilibrium distribution is realized only in the asymptotic 
limit t — )■ +00, in which there is no dependence on time and on the initial 
data. 

The potential V{(p) from ffTOj) can be written as V{ip) = 6S{(p)/6{p, where 
the static action functional is 

^oT/ f (v^W - </'(y))^ 



RdxRd Rd 

(14) 

in the d - dimensional coordinate space with coupling constant Qq. For proof 
of the statement V{ip) = 6S{ip)/6{p it is enough to write the difference S{(p + 
5ip) — S{ip) and to take into account of the symmetry with respect to values 
X and y in f lT4|) . It's important, that the first term in f lT^ is nonlocal, i.e. 
which is not polynom or function of the field or its derivative evaluated at a 
single point in the space. This reflects physically the existence of a long-range 
interaction in the system, and confirms the choice of operator generalization 
of the spatial derivative ([6]). Near a critical point all interactions are long 
range, it is confirmed by an experimental data. 

Application of the renormalization group method proposes the transform- 
ing of some parameters on renormalized. The non-renormalized parameters 
are marked by the subscript "0" to distinguish it from corresponding renor- 
malized. The functional f[T^ at m = 3 is nonlocal generalization of the 
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Ginzburg-Landau functional, in the field theory terms it is the massive 
model with mass 

For the renormalization of Eq. f fTOj) we will use the theorem of equivalence 
of stochastic dynamic and quantum field theory with two-component field 



with the main field y9(x, t) and the auxiliary field (p'{'x,t) 19|]. In the proof 
of this theorem assumption about just first-order of time derivative was not 
used. It was required only linearity of time-dependent operator. Therefore 
Eq. (ITU!) is equivalent to field - theoretic model with action 

+00 

S{<^) = j rf^x j dt{\^^'^' + ^'{-lDtv + ^oVm) (15) 

Rd -00 



where $ = Lp'). For the Cauchy problem ([8]) the action functional is 

+00 

5($) = y d'^x j dt{Xo^'^' + {-lD^^ + XoV{if) + 5{iL)5{t-t'))) (16) 

Rd 

that is reduced to adding to the action of term (y9'(x = 0, t = 0) , since always 
t > t'. Statistical correlation functions (Green functions) are defined by 
functional average with weigh exp S'($) (we include in S the usual minus sign 
in the argument ). Further we will work directly with the field - theoretical 
model fllSp . We are concerned with the restriction of parameters a,a,d, 
allowing to apply the renormalization group method. 

The asymptotics of correlators and other quantities at r — > +0 with al- 
ready the renormalized parameter r is called the critical asymptotics I?!. The 
physical meaning of this parameter is the deviation of the critical temper- 
ature: T T — Tc- Our model ata— flja— )'2, m— j-S transforms to 
the model A, where unique non-trivial critical asymptotic is dependence of 
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relaxation time on r. It is determined by the critical exponent z and has 
the form r^^^ . Other critical asymptotics coincide with those in the static 
model. For arbitrary a, o and m it is false. In following section the dif- 
ference between statics and dynamics for the field - theoretical model f lTSj) 
with fractional derivatives for arbitrary vertex degree m will be shown. 

3. Determination of propagators 

Standard Feynman diagram method with vertex (/'V^ bare propa- 
gators is corresponded to action fllSp . These propagators are found by a 
representation of the quadratic (called "free") term of the action 

+00 

5o($) = jd"^ j dt (Ao</^V + ^' - Ao(-A)'^/V - Aoro</^)) (17) 

Rd -oo 

as the integral of a quadratic form with a symmetric operator matrix K G 
M(2,2): 



hoo 



^o($) = -^j d'^^ J dt^K^ (18) 

Rd 00 

where $ = {ip,ip'). Then K^^{'x — x.',t — t') is the translationally invariant 
matrix of functions such that KK~^ = /5(x — x')5(t — t'), where / is the 
identity matrix. Using integration by parts rule for fractional derivative with 
respect to t, we obtain from f lT7|) the following matrix expression of matrix 
K in the momentum representation (let S = Xo{k'^ + tq)): 



'ID- + E 
IDf + E 2Ao 



(19) 



11 



where ^D^ is right Caputo derivative defined as 




n 



(20) 



We represent K ^ as 



^n(x-x',t-t') ^i2(x-x',t-t') 
^2i(x-x',t-t') 



) 



(21) 



where the functions Qu and Q12 called propagators are sought in the class 
of functions vanishing at ||x — x'|| — )■ 00. For correctness of the iteration 
procedure for the equation ( ITOi) and respectively the diagram technique it 
is also necessary, that the condition of stability Q12 — )• at t — t' — +00 
is complied. Typically, the propagator is normalized to delta - function at 
coinciding times, and the requirement of delay corresponding to the physical 
principle of causality is imposed [3, Q • It can be expressed by 




From this we see that to find the propagator Q12 we can use the Laplace 
transform in the variable t — t' to the variable s and the Fourier transform in 
the variable x — x' to the variable k. For convenience, different representa- 
tions of propagators differ in this text by their arguments. Then the equation 
KK~^ = /5(x - x')S{t - t') with ([19]) and ([2l]) - ([22]) can be rewritten as 
follows at t > 0: 



We use the property of the Mittag-Leffler function Ea{ut") for Re{s) > 





(22) 



s"6?i2(k, s) - s''-'gi2{k, +0) + Eg,2ik, s) 



(23) 



u 



l/a. 



12 



/ e'''Ea{ue)dt 



.a-l 



U 



(24) 



We finally obtain from ([22]) - (EH): 



Qi2{t - 1', k) = e{t - t')E^ {-E{t - t'Y) 



(25) 



Similarly, we obtain the advanced propagator: 



g^^it - 1', k) = e{t' - t)E^ {-E{t' - tT) 



(26) 



Note, that the propagator (125|) coincides with solution of the Cauchy problem 
(IZl), supplemented by the condition t > 0. Just as in model A, Q12 is the 
retarded propagator and Q21 is advanced. The propagator Qu can be uniquely 
determined by the convolution of functions Q12 and ^21- From ([T9|) .( l2T|) we 
obtain: 
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+00 



gn(K t) = 2Ao j e?i2(k, t - t')g2i(K t')dt = ( 

— oo 

+00 

= 2Ao J e{t - t')e{-t')E^{-E{t - t'r)Ea{-^{-tT)dt = (: 
—00 

(t > 

e{-t) J E^{-E{t-tT)Ea{-^{-tT)dt + e{t) J E^{-E{t-tT)Ea{-^{-tT)dt 
—00 —00 y 

+00 

= 2Ao J E^{-Et"')E^{-E{t' + \t\r)dt^i 


+00 

= 2Ao5-^/" J E^{-t"')E^{-{t' + \t\E}'''Y)dt{^ 


Usually the propagators Qu and Q12 are derived in the frequency - mo- 
mentum representation. If a is fractional, the time derivative is replaced by 
factor (ia;)", this factor fixes a single-valued branch. This is a well-known 
problem of non-physical poles separation. In our view it is more convenient 
to use the (k, t) - representation. 

4. Renormalization of field theoretical model in one-loop approxi- 
mation 

4-1. The analysis of scale invariance 

We consider in this section the more general class of models with power- 
law nonlinearity c/?"^: 

+00 

S{^) = - fd"^ f dt^^^ (32) 

J J ml 

Rd -00 
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It is convenient to do stretching variables in the auxihary field ^p' — )■ (/j'cq 





with Co = Aq ^. Than the action (l32l) takes the form 



+00 

m 



+00 

Rd -00 



(33) 

We will consider the theory in the coordinate space of dimension d .Fur- 
ther it is necessary to analyse scaling dimensions. We introduce the definition 
of two independent canonical dimensions: the momentum dimension /x^, the 
'requency dimension [i^ and total for an arbitrary value a [la = a/x^ + a^i^ 
3. By definition, /i^ = -/i^ = /i^ = -^^ = 1, /i£ = /i? = /i^^ = /i^ = 0, 
and dimensions of other quantities in f|T5|) are found requiring dimensionless 
of each term in the action. From this follows, that scaling dimension of in- 
tegration over all space and time is equal a + ad, scaling dimension of time 
derivative is a, the scaling dimension of spatial derivative is a. Then the 
system of equations, expressed dimensionless of the action ( !T5|) is 

/ico + = ad + cr (34) 
A^co + A*,^ + yU<^' + ao" = ad + (J 
[i-^p + + OLG = ad + a 
f^'fi + /^To + fi,p' = ad + a 
fJ'ifi' + l^go + T^ii^ = ad + a (35) 

Considering jointly the second and third equations in (1551) we get fico = 0, 
from third and fourth we obtain firo = ol^- We obtain immediately from first 
equation in fl35l) dimension of the auxiliary field /i<^/ = [ad^o^jl. Dimension 
of the main field is obviously [i^p = {ad+o')/2 — aa. Finally, the fifth equation 
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follows 2fig = ad{l — m) + cr(l — m) + 2aam. 

Further, the theory is considered in the space, where the coupling constant 
go becomes dimensionless. In this case existing divergences in correlators, if 
the condition of multiplicative renormalizability is complied, could be man- 
ifest in the form of poles at the deviation of the space dimension e. This 
space dimension d* is called logarithmic. In our case 



Despite the possibility of analytic continuation of diagrams in the region 
d < 0, we think, that it is not correspond to the physical meaning of the 
problem. Equation (!36ll clarifies restrictions on parameter a as condition on 
positiveness of spatial dimension. Considering ( l36l) . they become 



This inequality (I57|) along with the condition of the fractional representation 
of the Laplacian as the principal value of the integral operator and (1131) 
1 < 0" < 2 is conditions for the applicability of the quantum-field renormal- 
ization group method to the equation (|T0|) . Hence we see as increasing degree 
of nonlinearity m in the action narrows the range of permissible a. The equa- 
tion ( 136|) expresses the relation between the dimension of the coordinate space 
d and the orders of derivatives a and (3. This provides a self-consistency of 
theoretical-field model (132 p . From inequality (m — l)/2m < 1/2 implies that 
for the usually considered range 1/2 < a < 1 always satisfies ( P7|) . 

Now for applying the renormalization group method in the space di- 
mension d = d* — e it is necessary to prove the multiplicative renormal- 
izability of the theory p^]. The form of required counterterms is found 





(36) 
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using scaling analysis of 1-irreducible Green functions. Such functions with 
Nis> = N^p, N^i external legs have scaling dimension 5 = ad + a — ^q,Nq, = 
ad + a — fi^p'N^i — fi^pN^, where is field scaling dimension. The divergence 
is determined by the total canonical dimension 6: the diagram has surface 
divergence, if in logarithmic theory 6 is integer non-negative number. In our 
case at e = we find from (!35|) . (!36|) 

S(m — 1) , , 

^ = 2m-N^- mN^, (38) 

aa 

This means that surface divergences exist only in 1-irreducible functions 
{(p'ip'), {^p'^p) and {{p'ip"^), i.e. only in the functions presented in f lTSj) . 
Now we can write the renormalized action functional 



+ 00 



gM'ip 

5" 



m! 
(39) 

with renormalization mass M, renormalized parameters c, r, g and relations 
between dimensionless (they depends on g only) renormalization constants 
Z: 

Zi = Z^piZc 

z, = z]l^z]i-z^ 
= z]l^z]i- 
z, = z'J'z'J'z^ 

Z, = Zl('z:^/'Z, (40) 

It is known, that in model A there is a single new renormalization con- 
stant Zc (compared with the static case). Equality of other constants follows 
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from uniformity of Schwinger equation for static and dynamic generating 
functionals. The proof of this is based on V{ip) = 6S{(p)/6ip, where S{(p) is 
the renormahzed static functional. 



S{ip) = -Z^-P.V. J \\^-y\\d+a d^dy-J l^Z^Z^-if + ZgZ'^ ^— ^ 

RdxRd Rd 

(41) 

In general case we have following 

Theorem 1. The equality of the renormalization constants (^ip and / T^gj) 

Zip = Zipi = Zip^ Zg = Zg, Zj- = Zj- (42) 
holds if and only if a = 1. 

The proof is constructed like model A, taking into account the fact that 
the logarithmic dimension of the static action functional ( |T4l) is 

,* m + 1 

m-l 

Obviously, the coincidence of the logarithmic dimension f HSj) and f l36|) takes 
place only at a = 1. Thus, all UV (at large momentum) divergences mani- 
fest themselves as poles at e in the renormalization constants Zj, i = 1...5. 
Further, standard renormalization group equations expressed the renormal- 
ization invariance are written. We introduce the solution of Eq. ( ITOj) as 
follows at Co = (to, Ao): 

G'(x,t|eo) = J Vip' j Vipip{x,t)exp{S{ip',ip)) (44) 

Here the normalization constant is included into the differential Vip'Vip. 
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m + 1)!^ 



The condition of multiplicative renormalizability implies the relation G'(eo) = 
Z^^Gni^e, M) for the corresponding Green functions in Eq. (jH]). We use T>m 
to denote the differential operation MSm for fixed cq and operate on both 
sides of this equation with it. This gives the basic RG equation V^Gr = 0: 

[Dm + (3{g)dg - -fg]GR{e, M) = 0, = xd, (45) 

where the anomalous dimension 7^ in using the minimal subtraction scheme 
(MS) is expressed in terms of the renormalization constant Zg-. 7g(5') = 
Vjif^^Zg = —egdgln Zg/{1 + gdglnZg) and doesn't include poles at e. We 
need to further analyse the expression for the Gell-Mann-Low function P{g) = 
T>^fg in terms of the anomalous dimension of the coupling constant g: 

f3{g) = -g{s + ^gig)) (46) 

The anomalous dimension of an arbitrary value a expressed in terms of the 
corresponding renormalization constant as 7^ = f3{g)dg\n Za- 

Obviously, that only as a = 1 and the logarithmic space dimension of 
the dynamic theory coincides with the static case, we can use theorem for 
the coincidence of all renormalization constants excluding Z^ and accept 
behaviour of the Gell-Mann-Low P{g) function of the static model p^. In 
the general case of the dynamic problem (fT5|) we can not do it. 

4-2. Calculation of the one-loop diagrams 

Now we consider the case m = 3. In one-loop approximation there is only 
one diagram v^V^ (Fig- [I]) ? which contribute the ordinary term in renormal- 
ization constants: Zg^Z^-. We can calculate it in the MS-scheme at the null 
external momentum and the null external frequency. The condition cjgit = 
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Shi 

Figure 1: The only one-loop diagram, propagators Qn and Q12 are signed. 



is equivalent to the equality of time arguments of propagators. The pole part 
of this diagram T^>^p3 is denoted as xl^^cr): T^^v^ ~ x(q^)'^)/^- 

In MS scheme only the function x{<^j ^t) interest us. For its calculation 
in diagram all external momentum are supposed zero. The diagram T^/^s will 
equal: 

00 00 

r^V^ = J^^^J dkk''*-'-^gnik,t)guik,t) (47) 


where the logarithmic dimension d* is defined by fl36l) . Also we omit factors 
independent on a and a in the diagram. The final factor is determined from 
a compliance with the static case. Set in ( H71) ( l25l) and ( l3Ti) . we can submit 
this as an integral: 

00 00 00 
r^,^3 = JdtJ dkk''*-'-'E^{-Er) J c/rE,(-r°)E,(-(r + tS^/")") (48) 
00 

Replacing Et" — )■ we then obtain the product of two integrals 
r^V^ = 2A / dk'^—— [dt [ dTE^{-r)E^{-T^)E^{-{t + rr) (49) 
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The integral over the momentum space is now 



oo 



dk^^ = r^^''-'''^--\-- Jdyy 







'T a ^ ' cy. \ a. 



-^(.-.-i)-i^-#r(f-3 + f)r(3-i-f) 



and the diagram ( HHIl is equal 



W = r-^^" -'-'^-^X-^^ f;^ ^^J(«) (51) 



1 2r(^-3 + ^)r(3-i-^ 

-(a — e— 1) — — \ — — ' (t' V a a 

"A a _ , n , 

where we denoted 

oo oo 

J(a) = y y rfrE^(-t°)^a(-r")Ea(-(t + r)°) (52) 



In particularly at a = 1 Ea{x) = exp(x) and J (a) = 1/4. 

The expression fl50l) has a simple pole in e at a = 1 and atQ; = l/r;,, riG 
Z, n ^ 3. However, the situation n ^ 3 is forbidden by the condition 
( 137|) . Therefore at a = 1 and 1 < o" < 2 we have x(1,(t) = (t. In a range 
1/3 < « < 1 the expression (l50l) is regular, thus x(«, o") = 0. It is convenient 
to construct the perturbation series of the coupling constant g' = g/Sd, 
where 5"^ - surface area of the d - dimensional sphere: Sd = /2'^T{d/2) 
[3]. Than in the case < a < 1/3 the renormalization constant Zg will be 
regular in the one-loop approximation. Hence, if 



%{9')=lf^9' + 0{g'') (53) 
then 7^^-' = 0. But in the case a = 1 we conclude 

7« = Ca (54) 
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where C is a constant which is not depend on a. It is known, that in the 
model A at ci = 4 P^{g') = -g'e + 3g'^. Finally from g6]) and (ES]) , ([Ml) we 
conclude 

(3{g') = -g'e+^ag'' (55) 

The important consequence of the formula fl55l) is the existence of the IR 
stable point of the Gell-Man-Low function g'* = 2e/3a ~ e : (3{g'*) = 0. 
This fact confirms the critical scaling ata = l, 0<cr<2in generalized 
model A. 

5. Conclusion 

We applied field theoretic renormalization group tool to the stochastic 
fractional derivative nonlinear equation f llOp . In principle, we could real- 
ize this procedure even for the nonstochastic Cauchy problem but it 
is important for us to demonstrate a possibility of the MSR formalism ap- 
plication to fractional processes. Our work represents a generalization of a 
critical dynamics research observed in the model A. Inclusion of fractional- 
order derivatives of a with respect to time and of cr with respect coordinates 
to the model leads to some difficulties. 

Replacing the standard Laplacian by the fractional analogue leads to a 
nonlocal both action functionals: static and dynamic. That is why the prop- 
agators are sought in the momentum representation, in which the action is 
local. There are not other difficulties due to the substitution A — > — (— A)'^/^ 
apart a nonlocality of action, and the generalization is trivial. The case 
a = 1, cr < 2 is called super diffusion. 
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Propagators of our model at a 7^ 1 are expressed in terms of the Mittag- 
LefHer function Ea{x) depending on the complex parameter a. The propa- 
gator fl3Tl) has a complex integral representation for multi-loop calculations. 
This is due to the fact that the semigroup property of the Mittag - Leffler 
Ea{—x"') + Ea{—y°') = Ea{—{x + y)"') takes place only at a = 1. In this case 
the Mittag - Leffler function coincides with the exponential function. 

We prove the multiplicative renormalizability of the model with the non- 
linearity of the general form if"^ and found its logarithmic dimension (136 p . 
In the physically interesting range 1/2 < a < 1 the logarithmic dimension is 
always positive. 

One application of the renormalization group method to the theory of 
phase transitions is the establishment of scaling behaviour in the infrared 
asymptotic region /c — )■ 0, a; — )■ 0. A necessary condition for this is the 
fact that the Gell-Mann-Low function has a fixed point: f3{g*) =0, (7* ~ e. 
From the one-loop calculation we can see, that such situation is realised with 
the restrictions fl37|) only for the superdiffusion a = 1, a < 2. For the 
subdiffusion a < 1, cr = 2 we can not assert something specific. 

Accordingly if a = 1 and the spatial dimension is less than d* = 2a, 
then in the systems controlled by the equation f lTUlIT^ the critical scaling is 
predicted by us. 

6. Acknowledgements 

This research is supported by the Chebyshev Laboratory (Department 
of Mathematics and Mechanics, St. Petersburg State University) under RF 
Government grant 11. 034.31. 0026 



23 



References 

[1] M. B. Isichenko, Rev. Mod. Phys. 64, 961 (1992) 

[2] G. Falkovic and K. R. Sreenivsan, Phys. Today 59(4), 43 (2006) 

[3] V. E. Tarasov, G. M. Zaslavsky, Physica A, 368:2 (2006), 399-415 

[4] Thermodynamics: An Engineering Approach ((Fourth Edition) ed.). 
McGraw-Hill. 2002. pp. 824. 

[5] De Dominicis C., Brezin E., Zinn- Justin J., Phys. Rev., B12:ll (1975), 
494-495 

[6] Hohenberg P. C., Halperin B. I., Rev. Mod. Phys., 49:3 (1977), 435-479 

[7] A.N. Vasiliev, Functional Methods in Quantum Field Theory and Sta- 
tistical Physics, Gordon and Breach, New York, 1998 

[8] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, Oxford 
Univ. Press, New York, 1989 

[9] E. Lutz, Fractional Langevin equation, Phys. Rev. E 64 (2001) 051106- 
1-051106-4 

[10] O. P. Agrawal, Solution for a fractional diffusion- wave equation defined 
in a bounded domain. Nonlinear Dynamics 29 (2002), 145-155 

[11] Antonov, N. V, Field-theoretic renormalization group for a nonlinear 
diffusion equation. Physical Review E, 66:4 (2002), 046105 



24 



[12] E. Barkai, Fractional Fokker-Planck equation, solution, and application, 
Phys. Rev. E 63 (2001) 

[13] Renormalization of tracer turbulence leading to fractional differential 
equations, Phys. Rev. E 74 (2006) 

[14] N. V. Antonov, A. N. Vasiliev, Critical dynamics as a field theory, Teo- 
ret. Mat. Fiz., 60:1 (1984), 59-71 

[15] J Bonart, L F Cugliandolo and A Gambassi, Journal of statistical 
mechanics-theory and experiment, P01014, 2012, Critical Langevin dy- 
namics of the 0(N)-Ginzburg-Landau model with correlated noise 

[16] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and Apphcation 
of Fractional Differential Equations, North-Holland Math. Stud., 204, 
Elsevier, Amsterdam, 2006 

[17] Gorenflo, R. and Mainardi, F. (1998), Random walk models for space- 
fractional diffusion processes. Fractional Calculus and Applied Analysis, 
1, 1677-191. 

[18] Hanyga, A. (2001), Multi-dimensional solutions of time-fractional 
diffusion-wave equations, Proc. R. Soc. London A. 

[19] Martin P. C, Siggia E. D., Rose H. A., Phys. Rev., A8:l (1973), 423 - 
437 



25 



